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ABSTRACT 



Aims. A long-timeline kinematic study of the archetypal compact symmetric object (CSO) OQ 208 sheds light on the 
physical properties of the most compact radio sources. 

Methods. Archival data from the Very Long Baseline Array (VLBA) at 15 GHz over a time span of 13.6 yr were used 
to investigate the kinematics of the radio source. The flux density monitoring data obtained at the Michigan 26-meter 
radio telescope were also used as supplementary information for analyzing the geometry of the radio structure. 
Results. At 2.3-GHz, the radio emission is dominated by two mini-lobes separated by ~10 pc in a northeast-southwest 
(NE-SW) direction. At 8.4 and 15 GHz, each lobe is further resolved into two subcomponents, which are identified as 
hotspots. A knotty jet is linked with the NE hotspot and traces back toward the geometric center. The core is too weak 
to be detected. Significant flux density variation is found in the primary hotspots with a maximum level of 62% (NE1) 
and 19% (SW1). The flare epoch of NE1 is earlier than that of SW1 by approximately 5.00 yr, suggesting that the 
northeast lobe is advancing and the southwest lobe is receding. This light travel difference indicates a radial distance 
difference between the two hotspots of 1.53 pc, which indicates an inclination angle of about 80.8 degrees between the 
radio jet and the line of sight. The angular separation rate between NE1 and SW1 is 0.027 mas yr -1 , corresponding to 
a projected speed of 0.133 c. The inner jet knot (Jl) moves at 0.047 mas yr~ (or 0.230 c), about 3.5 times the hotspot 
advancing speed. 

Conclusions. The large viewing angle and the modest jet speed suggest a mildly relativistic jet. The jet axis is close to 
the plane of the sky. The separation rate and the distance between the two primary hotspots result in a kinematic age 
of 255 ± 17 yr, confirming that OQ 208 is indeed a young radio source. In addition to the hotspot advancing motions, 
sideways motions provide evidence that the lobes are obstructed by the external interstellar medium. 

Key words, radio continuum: galaxies - galaxies: active - galaxies: individual: OQ 208 



1. Introduction 

Compact symmetric objects (CSOs) are a subclass 

of extragalactic radio sources that are characterized by 
a compact double or triple radio structure with an 
overall size less than 1 kpc. The physical nature of the 



are confined within the host galaxy and experience stag- 
nated growth. The CSO ages provide a critical distinction 
between the two models. Previous kinematics studies of 
individual CSOs and subsamples of these sources (e 



youth model (e.g. |Phillips et al. (1982) 


Fanti et al. (1995)1 


Readhead et al.(1996)p proposes tha 


t CSOs are small 



Owsianik & Conway (1998)1 


Taylor et al. (2000) 


qow that the 


An et al. (2012)} s 



Owsian ik et al. (1 998 
[Gugliucci et al. (2005' 



tic radio source evoluti o n; the frustration model (e.g 



van Breugel et al. (1984)} |Q'Dea et al. (1991)[ ) attributes 
the small size of CSOs to extremely strong confine- 
ment by the dense external medium. The two models 
define distinctly different evolutionary fates of CSOs. 
In the youth model, all young CSOs eventually evolve 
into large-scale double sources, i.e., Fanaroff- Riley (FR, 



tter from ~ 0.04 c to 
~ 0.5 c, resulting in young kinematic ages in the range of 
only 100-2000 yr. However, additional sideways motions of 
hotspots and disturbed lobes indicate strong interactions 
between th e jet heads and the sur ro unding interstella r 



medium (Stanghellini et al. (2009) 
Sideways motion is usually very sio~ 



An et al. (2012)[ ). 
compared to the 



Fanaroff & Riley (1974) ) type-II sources, over a few million 
years; whereas according to the frustration model, CSOs 



dominant advancing motion, which requires highly 
accurate measurements. A detailed investigation of well- 
selected CSO samples is essential for exploring the complex 
kinematics of CSOs, and for understanding the physical 
environment of the host galaxies on 1-500 pc scales. 
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The kinematic age (r^) of CSOs is traditionally deter- 
mined by dividing the separation (R) between two termi- 
nal hotspots by the separation rate (/i) . A high-accuracy 
measurement of requires high-precision position deter- 
mination for individual observations, and a sufficiently long 
timespan of the data. For more than two epochs, /i is of- 
ten determined from a linear regression fit to the chang- 
ing hotspot separation with time. The statistical uncer- 
tainty of the fit is sensitive to the number of available data 
points and the uniformity of the time sampling. In addition 
to this, different resolutions of the interferometric images 
and the intrinsic opacity effect in different levels may in- 
troduce systematic errors in the proper motion measure- 
ments. An accurate measurement for the separation rate 
/x requires CSOs with multiple-epoch Very Long Baseline 
Interferometry (VLBI) imaging data at the same observing 
frequency over a sufficiently long time span. 

In this paper, we carry out a kinematic study of an 
archetypal CSO, OQ 208, on the basis of archival VLBI 
data over a time baseline of 13.6 yr. OQ 208 (also known 
as Mrk 668, B1404+286) is one of the closest CSOs (z = 
0.0766: Huchra et al. (1990) I and provides a template for 
the dynamic properties of the most compact (and possibly 
the youngest) radio sources. The host galaxy of OQ 208 
shows typical Seyfert-1 spectra with strong broad Balmer 
lines (FWHM ffo = 6000 k m s" 1 ) and also forbidden lines of 
[Ne III], [O III] and [S II] dBurbidge fc Strittmatter (1972) 



Eracleous &; Halpern (1994)p . In the radio band, a broad 
21-cm HI absorption line was detected in OQ 208, 
indicating a fast and massive outflow of neutral gas 
dMorganti et al. 2005[ ). The radio emission is concen- 
trated in a compact ce ntral region within ^10 mas 
(Stanghellini et al. (1997) I. Additional evidence for the 
compactness comes from the convex radio spectrum with 
a turnover at about 4.9 GHz (Dallacasa et al. (2000)), re- 
sulting fr om synchrotr on self-absorption or free-free ab- 
sorption ( |Fanti (2009)] ) . The VLBI images of OQ 208 
at 2.3 and" 5 GHz reveal a typical CSO-type mor- 
phology with two compact (mini-) lobes along the NE- 
SW direction ^Stanghellini et al. (1997)[ |Liu et al (2000)| 
Wang et al. (2003)p . At higher frequencies of 8 and 15 GHz , 
a more detailed structure is revealed: the mini-lobes are 
resolved into subcomponents, a nd internal jet knots are 
detected between the two lobes ( Kellermann et al. (1998)] 
Luo et al. (2007)1 Previous proper motion measurements 



of OQ 208 (e.g-ILiu et al "2 000) , Stangh ellini et al (2002) 



Wang et al. (20037]" Luo et al. (2007)P were mostly base< 



on 5- and 8-GHz data with lower angular resolutions and 
are affected by the opacity effect. Moreover, these measure- 
ments only cover a time span of a few years, insufficient for 
tracing the long-term change of the hotspots. In the present 
study, we only made use of 15-GHz VLBI data to deter- 
mine the separation velocities of hotspots and the proper 
motions of the internal jet knots. 

The structure of the paper is as follows. Section [2] 
describes the VLBI data used for the kinematics study. 
Section [3] presents the results, including radio images, light 
curves, and proper motion measurements. A conclusion and 
a summary are given in Section [4j Throughout this paper, 
we assume a flat cosmological model with Hq = 73 km s _1 
Mpc" 1 , il M = 0.27, n A = 0.73. At the redshift of OQ 208, 
1 mas angular size corresponds to 1.4 pc projected linear 
size, and a proper motion of 1 mas yr -1 corresponds to 4.9 
c apparent speed. 



2. VLBI data 

Owing to its compactness and high brightness, OQ 208 
is often used as a calibrator in VLBI experiments. 
Consequently, there are abundant archival VLBI data sets 
from the past two decades. The 15-GHz VLBA data from 
the Monitoring of Jets in Active Galaxies with VLBA 
Experiments (MO JAVE) program^ have the highest resolu- 
tion and sensitivity and are most appropriate for kinematic 
studies. The 15-GHz data spread over 28 epochs in the 
time range from early 1995 to the middle of 2009. We also 
made use of the VLBI data from the Radio Reference Frame 
Image Database (RRFID)0, which were obtained with the 
VLBA together with several geodetic radio telescopes si- 
multaneously at 2.3 and 8.4 GHz. The multiple-frequency 
data were used for the spectral index analysis of compact 
components. The RRFID data cover the time range from 
1994 to 2008. 

The calibration of the archival VLBI data has al- 
ready been made in the Astronomical Image Processing 
System (AIPS) following the standard procedure. We per- 
formed additi onally several iterat ions of self-calibration 
in DIFMAP ( |Shepherd et al. (1994)] ) to eliminate resid- 
ual antenna-based phase errors. The 15-GHz visibilities 
were fitted with five Gaussian model components (two in 
the northwest lobe, another two in the southeast lobe, 
and one jet knot, see Figure [1} using the DIFMAP task 
MODELFIT. Table 1 lists the fitted parameters, includ- 
ing the integrated flux density Sj, the relative separation R 
from the primary hotspot used as the reference, the posi- 
tion angle of the VLBI component, and the deconvolved size 
6maj x dmin- The statistical uncertainties of the observed 
pa rameters were esti mated according to the formulae given 
by Fomalont (1999)| except that we also included an addi- 
tional 5% as the amplitude calibration error. The details of 
the error analysis method are given by | An et al. (2012)] 

In addition to VLBI imaging data, the single-dish 
flux density monitoring data observed with the 26-meter 
paraboloid telescope of the University of Michig an Radio 
Astronomical ObservatorjQ (Aller et al. (1985)1 were in- 
cluded for a supplementary analysis of the variability and 
the geometry of the radio structure. The total flux density 
measurements were made at 4.8, 8, and 14.5 GHz from July 
1974 to August 2009. 

3. Results 

3.1. Radio morphology 

At 2.3 GHz, OQ 208 exhibits the double-component mor- 
phology typical of the CSOs (Figure QJa). The total extent 
of the radio source is about 7 mas (~10 pc). The northeast 
(NE) component is much brighter than the southwest (SW) 
one with an intensity ratio of 18.0:1 at 2.3 GHz. Both com- 
ponents show a steep spectral index at high frequencies 
(defined as S oc v~ a ) with afj^ z = 1.38 ± 0.03 (SW) 
and al 5 gf H z z = 1.12 ± 0.02 (NE), identifying them as two 
mini-lobes of the CSO. The asymmetric brightness of two- 
sided jets/lobes in radio sources is commonly attributed to 
the Doppler boosting effect that enhances the apparent flux 
density of the advancing jet/lobe by a factor of <5 3+Q , where 



1 http:/ /www. physics. purdue.edu/MOJAVE/ 

2 http:/ /rorf. usno.navy.mil/rrfid.shtml 

3 https:/ /dept. astro. lsa.umich.edu/datasets/umrao.php 
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Fig. 1. Total intensity image of OQ 208. Panel a shows the 2.3-GHz VLBI image, displaying a double-component 
morphology. Panel b shows the tapered 2.3-GHz image. The extended emission at ~30 mas west of the primary structure 
is likely a relic radio emission resulting from the past activity a few thousand years ago. Panels c and d show the 8.4 
and 15 GHz images. Each mini-lobe is resolved into two hotspots. Between the hotspots NE1 and SW1, a knotty jet is 
detected. 



5 is the Dopplcr boosting factor and a is the spectral in- 
dex. However, the kinematic analysis below suggests that 
the jet in OQ 208 is mildly relativistic and the jets nearly 
align with the plane of the sky, therefore Doppler boosting 
cannot account for the brightness difference. Another expla- 
nation of the strong asymmetry in the brightness of the 
two lobes is that the receding SW lobe suffers from more 
free-free absorption than the advancing NE lobe. Indeed, 
the flux density ratio is even hi gher at a lower freque ncy of 
1.66 GHz, S N e/Ssw = 60 : 1 ( |Kameno et al. (2000)[ ), sup- 
porting this interpretation. A third possibility is that the 
inhomogeneous distribution of the external medium sur- 
rounding the lobes results in a larger conversion efficiency 



fro m jet kinetic energy to radiat ive energy in the northeast 
jet (|Orienti k, Dallacasa (2012)[). 



A 2.3-GHz image with a Gaussian taper with a half- 
value at 20 MA wavelength (FigurcQJb) reveals an extended 
feature about 30 mas ( ~40 pc) to the w est, which was pre- 

. Extended emission 
so been detected in 



viously reported by Luo et al. (2007) 
features on kpc to Mpc scales has a 



OQ 208 ( Ide Bruyn (1990)1) and i n other CSO galaxies (e.g., 
0108+388: |Baum et al. (1990)1 |Stanghellini et al. (2005)| 
0941-080, 1345+125: |Stangnellini et al. (2005)p , which were 
interpreted as relics remaining from past (> 10 8 yr ago) 
nuclear activity. Extended components on scales of <100 
pc are rarely seen in CSOs (J151 1+0518 is another ex- 
ample; Orienti & Dallacasa (2008) ) and this puzzling one- 
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sided extended feature at 40 pc distance requires an inter- 
val between two intermittent activities shorter than 2 x 10 3 



yr (Orienti & Dallacasa (2008)). The non-detection of the 
northeast fading lobe likely indicates asymmetric properties 
of the ambient interstellar medium (ISM) on pc scales, 
leading to more rapid radiative or adiabatic losses and a 
shorter life of the NE lobe. This, again, is consistent with 
the fact that the NE advancing lobe is much brighter. 

The NE lobe is resolved into two subcomponents (NE1 
and NE2) at 8 and 15 GHz. NE1 dominates the flux density 
of the whole source. The spectral index of NE1, determined 
from 8 and 15 GHz data at close epochs, is ctl^PJ 1 / — 0.99± 
0.18. The brightness temperature is calculated using the 
equation 



T b = 1.22 x 10 



12 



Sob 



ob nla J nlv 



-(1 + z), 



(1) 



where S b is the observed flux density in Jy, v Q b is the ob- 
serving frequency in GHz, 6 ma j and 9 m i n are the major 
and minor axis of the Gaussian model component in units 
of mas, and Tb is the derived brightness temperature in 
source rest frame in units of Kelvin. The average bright- 
ness temperature of NE1 is 4.4 x 10 10 K. The secondary 
component NE2 is weaker than NE1 in the range 4.3-16.3. 
It has a lower brightness temperature of 1.4 x 10 9 K and 
a much steeper spectral index of a^^f/ = 2.09 ± 0.13. 
The high brightness temperature and the relatively flatter 
spectral index of NE1 identify it as the primary hotspot 
formed by the reverse shock when the jet head impacts on 
the wall of the external medium. The continuous emission 
structure between NE1 and NE2 would suggest a physical 
connection. 

The SW lobe is resolved into two components (SW1 and 
SW2) at 8 and 15 GHz. The two components have compa- 
rable sizes. The integrated flux density of SW1 is slightly 
higher than SW2 with a flux density ratio R(SW1/SW2) 
in the range of 1.3-2.4. Their spectral indices are a^Qj^, — 
1.11 ±0.07 (SW1) and 1.82 ±0.26 (SW2). 

The presence of double hotspots in the 10-pc lobe 
structure of OQ 208 remains a puzzle. Similar mor- 
phology is found in the multiplc-hotspot appearance of 
large-scale double sources such as 3C 20 ( |Laing (1981)1 
Hardcastle & Looney (2001)), where the hotspot usually 
identified as the primary hotspot is extremely compact 
and bright, and the other is more diffuse and shows var- 
ious structures. Three models may account for the double 
hotspots inside the 10-pc nuclear region of OQ208: 

(1) Blocking of the jet head by the external medium. In 
this scenario, the jet impacts on the wall of the surrounding 
medium and generates the hotspot there. As the jet head 
may slide along the wall, the hot spot position als o changes 
accordingly (dentist's drill model: Scheuer (1982)). The pri- 
mary hotspot is the currently active jet-ISM interaction, 
while the secondary hotspot represents the past primary 
hotspot, which is now fading away. 

(2) Jet precession. The present images do not allow 
us to distinguish w hether a single jet beam is precessing 



Cox et al. (1991) I or if there are intrinsically two beams 
Laing (1982)) due to the similar compactness and shape 



of two hotspots and the absence of emission between the 
hotspots and the nucleus. One possible driving mechanism 
of jet precession is associated with a binary black hole sys- 
tem where the black hole launching the jet has a relative 



motion with respect to the surrounding ISM and the mul- 
tiple hotspots reflect the different impact locations. Taking 
the brighter SW1 as the primary (current) and SW2 as 
the secondary (older) impact of the jet beam(s) on the ex- 
ternal medium, the older hotspot has stopped receiving di- 
rect energy supply from the nucleus, and synchrotron aging 
will result in a steeper radio spectrum at high frequencies. 
Indeed, SW2 has a steeper spectral index than SW1. If SW2 
is the old hotspot, it should lie behind the primary hotspot 
SW1, implying that there is a jet bending between SW2 
and SW1. 

(3) Redirected flows. Alternatively, the secondary 
hotspot represents either a r e-directed out flow originat- 
ing from the primary hotspot dLaing (1981) I or a deflected 
jet flow ( |Lonsdale fc Barthel (19 86) ) escaping from a weak 
point in the cocoon/lobe structure. Different from the 
model (1), the secondary hotspot in model (3) continues to 
receive energy supply from the primary hotspot. Therefore, 
it would last longer, would not show spectral steepening, 
and be much fainter and less compact than the primary 
hotspot. In the eastern lobe, NE2 has no direct connec- 
tion with the nucleus, but shows a collimated bridge that 
links it with NE1. NE2 is also larger than NE1 and has 
an amorphous shape. These morphological characters make 
NE2 more likely to be the result of a re-directed outflow 
from the primary hotspot NE1 than a directed flow from 
the nucleus. 

The continuous emission between two hotspots NE1 and 
SW1 at 2.3-GHz (Figure [TJa) is resolved into a knotty jet 
at 8.4 and 15 GHz (Figure[Tj-c and[T]-d). From east to west, 
these knots are labeled Jl, J2, J3, and J4. Jl is the bright- 
est jet knot located at the inlet of the northeast lobe and is 
detected at 15 and 8.4 GHz at all epochs. It has a 15 GHz to 
8.4 GHz spectral index between 0.81 and 1.45, indicative of 
a typical optically thin jet knot, and a brightness temper- 
ature of 3.4 x 10 9 K. J2 appears in the 8.4-GHz images at 
epochs from 2002.044 to 2002.945, where the VLBI images 
have the highest sensitivity and resolution. J2 only appears 
in the most sensitive 15-GHz image at epochs 1995.958, 
1996.378, 2003.318, 2008137, and 2009.558. It has a steep 
spectrum with ct^Q^' = 1-43 and a brightness tempera- 
ture Tb < 1.9 x 10 K. J3 is closest to the geometric center 
and only appears in the highest sensitivity 15-GHz images. 
The brightness temperature of J3 is < 3.0 x 10 7 K derived 
from the 15-GHz data. Lacking the spectral index informa- 
tion, the active nucleus or the innermost jet knot nature 
of J3 remains uncertain. Regardless of the identification 
of J3, the non-detection at 8.4 GHz would imply an ex- 
tremely high absorption opacity toward the nucleus. J4 is 
only marginally detected at 15 GHz on 2009.558. More sen- 
sitive observations are necessary to to reveal the physical 
nature of these components. 



3.2. Variability of the hotspot NE1 

The steep-s pectrum lobes dominate the total flux density of 
most CSOs.lFassnacht & Taylor (2001) monitored a sample 
of seven CSOs using the VLA at 8.5 GHz and found ex- 
tremely stable flux densities (rms variability less than 1%) 
over a period of eight months. Some exceptional CSOs show 
moderate-level variability on time scales of a few years, 
for example, OQ 208 (the present paper) and J1324+4048 
(An et al. (2012) |. Unlike with blazars with a jet luminosity 
amplified by Doppler boosting, the low- and moderate-level 
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Fig. 2. Variability of VLBI components of OQ 208. (a): Total flux density variability of OQ 208 observed with the 
UMRAO 26-meter radio telescope: 4.8 GHz (black), 8 GHz (red), 14.5 GHz (blue). The sum of 15-GHz flux densities 
of the VLBI components (see Table 1) is represented by green symbols indicating that ~15% is missed due to lack of 
short uv-spacings; (b): Temporal variations of the flux densities of VLBI components. The flux densities are derived from 
the model-fitting results listed in Table 1. Two vertical lines mark the peak epochs of NE1 (around 1998.313) and SW1 
(around 2003.313). To clearly show the variability behavior of SW1, the flux densities of SW1 have been multiplied by a 
factor of 10; (c): Variation of the 15-GHz flux density of the hotspot NE1 with the component size. The flux density of 
NE1 experienced a slow rise before epoch 1998.313, followed again by a sharp decline until epoch 2005.473 followed by a 
slow decline. The component size constantly increases during the whole period; (d): Temporal variation of the spectral 
index a^J 1 / of the NE lobe. The spectral indices are calculated from the 15-GHz MOJAVE data and 8-GHz RRFID 
data at close epochs. A sharp jump of the spectral index occurs around 2006, after which the spectral index steepens 
significantly. 



variability observed in lobe-dominated CSOs is probably re- 
lated to variation of the energy supply by the relativistic jet 
flow and the changes of the energy dissipation process in- 
cluding the adiabatic expansion losses and radiative losses. 

The light curves observed with the Michigan 26-meter 
radio telescope (UMRAO) at three frequencies (4.8, 8, and 
14.5 GHz) show two distinct flares in the past 25 years, 
around 1983 and 2000 (Figure O-a). The accurate peak 
epoch around 1983 cannot be clearly determined because of 
the broad gap in the sampling, while the second-peak epoch 
is not well constrained because the light curves only cover 
the declining part of the flare. Th e 5-GHz flux density vari 



ation of OQ 208 using the VLA Stang hellini et al (1997) 
shows a continuous decrease from 1980 to 1995. In addition, 
Waltman et al. (1991) observed a decrease of the 8.1-GHz 
flux density from about 2.7 Jy in mid-1983 to about 2.0 Jy 
in mid-1989 using the Green Bank Interferometer. Starting 
from 1989, the flux density has been constant around 1.9 



Jy. This evidence consistently suggests that OQ 208 has 
intrinsic variability on time scales of a few years. 

For completeness, the 15-GHz summed flux densities of 
the VLBI components were added to diagrams covering 
the time range between July 1995 and August 2009, which 
confirms the second flare peak between 1998 and 2000. The 
declining part of the VLBI flux density data is fully consis- 
tent with the total flux density data. The comparison of the 
VLBI and single-dish data also suggests that a significant 
percentage (> 85%) of the total flux density originates 
from the compact radio structure. 

To quantitatively evaluate the variability, we intro- 
duced a vari ability index V to illustrate the relative vari- 
ability scale (|Zhao et al. (2012)[): 



V 



Sn 



Sn 



Sn 



(2) 
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where S ma x, S m j„ are the maximum and minimum flux 
density of an individual flare. A value of V « rep- 
resents indistinguishable or undetectable variability, and 
V ~ 1 indicates an extreme variability The 15- GHz light 
curves of individual components (Figure [H-b) suggest a 
high-variability index V = 0.62 ± 0.12 for NE1 during the 
1998 flare and a mild V = 0.19 ± 0.04 for SW1 during the 
2003 flare. The standard deviation of the measured flux 
densities with respect to the mean value is adopted as the 
statistical error. The variability index of NE1 derived from 
VLBI data is consistent with that measured from the single- 
dish monitoring data (Figure[5]-a) > confirming that the total 
flux density variation is dominated by the brightest hotspot 
NE1. The other VLBI components do not show significant 
variability at 15 GHz. 

The variability at 8 and 14.5 GHz using the single-dish 
data (Figure[5]-a) indicates a maximum level of (22 ±4)% 
and (40 ± 8)%, respectively. The 5-GHz variability scale 
is lower than at 8 and 14.5 GHz, probably because of 
the increasing opacity at lower frequencies. The mini-lobes 
at 2.3 GHz based on the RRFID data did not show any 
prominent variability. At 1.4 GHz, the total flux density of 
OQ 208 increased from —755 mJy in 1979 to ~854 mJy in 
2001 (dc Bruyn, private communication), which is closely 
connected with the growth of the overall source size. 

The complex variability behavior of OQ 208 can be un- 
derstood separately in the low- and high-frequency regimes: 

(1) Up to 1 GHz, the emission mostly arises from the 
extended structure surrounding the compact jets and lobes, 
which are strongly absorbed. The flux density of an opaque 
source can be simply expressed as S v oc 7], x A, where S v is 
the flux density at an optically thick frequency v 1 Tf, is the 
brightness temperature at frequency v 1 and A is the surface 
area of the extended structure. The continuous increase of 
the flux density from 1.4 GHz to 335 MHz results from the 
(slow) growth of the overall source size as a result of the 
expansion of the radio jets and lobes. 

(2) Above 1.4 GHz, the synchrotron radiation from the 
extended structure drops abruptly and the flux density 
is dominated by the compact hotspot, the mini-lobe, and 
the jet components. Arguments can be made that radia- 
tive losses play a dominant role in the energy balance of 
young and compact radio source s in the presence of strong 
(milligauss; mG) magnetic fields ( Stanghellin i et al. (1997) 



Orienti et al. (2006)||Orienti fc Dallacasa (2012)[ ). The ad! 
abatic expansion would then only lead to a decrease in 
the optically thin section of the spectrum. However, the 
source does not show prominent variability at 2.3-GHz, 
which could mean that the source is still opaque at this 
frequency and that any low-amplitude flux density varia- 
tion is smeared out by the opacity effect. 

(3) At frequencies above 5 GHz, the spectrum becomes 
optically thin and the observed variability is governed at 
the hotspots, i.e., the flux density from the most active jet- 
ISM interaction interface varies with the balance between 
the feeding and acceleration of fresh relativistic electrons 
and the radiative losses. The variation of the 15- GHz flux 
density with the component size (Figure [2}c) shows and 
increase of the flux density increases with the increasing 
size before the flare peak at epoch 1998.313, indicating the 
injection of particles, energy, and magnetic fields into the 
lobes. During this stage, the hotspot is still optically thick. 
After 1998, the intermittent feeding through the jet de- 
creases and the hotspot continues to expand adiabatically 



and to become more optically thin. The strong synchrotron 
radiation results in a sharp decrease of the flux density from 
~900 mJy at 1998.313 to -400 mJy at 2005.473. Figure^d 
shows the variation of the spectral index a^oifl* with time, 
using the 15-GHz MOJAVE data and 8.4-GHz RRFID data 
at close epochs. After the clear division around epoch 2006, 
the spectrum is significantly steepened. This supports the 
idea described above that the flux density variation of 
hotspot NE1 is regulated by intermittent nuclear feeding. 
When there are no or less fresh relativistic electrons in- 
serted into the hotspots, the old electrons suffer from syn- 
chrotron aging, and the spectral index in the optically-thin 
part steepens. 

Stanghellini et al. (1997) estimated the lifetime of the 
electrons radiating at 5 GHz assuming cquipartition mag- 
netic fields and a homogeneous spherical geometry for the 
hotspot. The resulting magnetic field in hotspot NE1 (com- 
ponent A in Stanghellini et al. 1997) is 90 mG, and an elec- 
tron lifetime at 5 GHz of 25 yr. The synchrotron aging time 
scales with the observing frequency and the brightness tem- 



perature Tf, as t syn 



2 T^. Using the same arguments, 



the lifetime of electrons radiating at 15 GHz is about 4 yr, 
in agreement with the observed variability time scale of ~15 
yr (from 1983 to 1998). 

3.3. Geometry of the radio structure 

According to Figure (5}b, the light curve of the northeast 
hotspot NE1 peaks around 1998.313. The peak of the light 
curve of SW1 is not very prominent and falls in a broad 
time range between 2002 and 2005. A mediate epoch of 
2003.31+t ;53 is chosen as an estimate of the peak epoch, in 
which the uncertainty is represented by the scattering of 
the peak epoch. If the radio flux density enhancement in 
both hotspots is associated with the same episodic nuclear 
activity, the observation that the flare of NE1 is lead- 
ing that of SW1 by b.OOtHl yr suggests that NE1 is 
the advancing lobe and SW1 is the receding one. The light 
travel time along the line of sight between the two hotspots 
corresponds to a radial difference of 5.00^53 light years, 
or 1.53j]Q'4y P c - Assuming that the two lobes are symmetric 
relative to the core (separation —9.5 pc) and that both jets 
are ejected at the same velocity, an inclination angle of the 

+2 8° 

radio structure of 80.8° " may be derived between the 

— 3.8 

jet and the line of sight. This calculation suggests that the 
jet axis of OQ 208 is very close to the plane of the sky. It is 
supported by the observation that the jet emission at the 
nucleus is not Doppler boosted (Sections 3.1 and 3.4). 

However, this implication appears to be conflict with 
the optical spectroscopic identification of OQ 208 as 
a Seyfert 1 galaxy with broad Balmer emission lines 



( Burbidge fc Strittmatter (1972) I. In the sta ndard unifi- 
cation model for AGNs ( |Antonucci (1993) I , the Type- 
I AGNs should have an inclinat ion angle less than 
70° ( |Zhang (2005)| |Liu et al. (2007)] ) . Such an angle of 
70° would predict a time delay between NE1 and SW1 of 
about 10.97 yr, for which there is no evidence in the light 
curves. The inconsistency between these inclination angles 
may suggest that the radio jet does not align with the op- 
tical axis of symmetry. 

If the brightness asymmetry of the two lobes is caused 
by differential free-free absorption (Section 3.1), the phys- 
ical properties of the ISM can be constrained. The spec- 
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Fig. 3. Relative proper motions of VLBI components with 
respect to the brightest hotspot NE1. The proper motions 
of the VLBI components are determined by linear fits to the 
separations as a function of time. The angular separation 
rate is calculated as p(Jl) = —0.034 ± 0.001 mas yr" 1 , 
H(SWI) = 0.027 ±0.002 masyr" 1 , and p(SW2) = 0.022 ± 
0.002 mas yr -1 . The hotspot NE2 does not show significant 
proper motion. 



tral fit of fre e-free absorption give s a differential opacity 
At// = 5.3 (Kameno et al. (2000)). The electron density 
is about 2700 cm -3 along a 1.53 pc path, assuming that 



the ISM surrounding the northeast and southwest lobes is 
homogeneous and the electron temperature is T e = 10 4 K. 
This electron density is typical for the narrow-line regions 
(NLR). An inhomogeneous ISM would require a higher den- 
sity (in a clumpy medium) to account for the free-free ab- 
sorption discrepancy. 



3.4. Kinematics and jet properties - a young radio source 

Existing kinematics studies of CSOs show that the hotspot 
advancing speed is typically around 0.2c, and that they are 
youn g with kinematic ages of on ly between 100 and 2000 
yr flOwsianik k Conway (1998)1 |Owsianik et al. (1998) 1 
Polatidis k Conway (2003)1 |Taylor et al. (2000)[ 



Gug iiucci et al. (2005)1 |An et al. (2012)p . In the present 
work, the slow proper motion of OQ 208 components are 
determined from making use of the long time-baseline 
15-GHz VLBA data to eliminate systematic errors induced 
by different resolutions, different uv coverage, and different 
frequencies. Since the identification of the central core of 
OQ 208 is not certain yet, the brightest hotspot NE1 is 



used as the reference, and relative proper motions of other 
VLBI components are measured. 

Figure [3] shows the variation of the distance of the VLBI 
components using 28 individual measurements over a time 
span is 13.6 yr. Linear regression was used to determine the 
separation rate. The positional uncertainty (a ~ 0.02 has) 
of each component at each epoch was used for weighting the 
individual data points (1/tr 2 ). The fits give the separation 
rates 0.027 ± 0.002 mas yr" 1 (0.134 ± 0.009 c, SW1-NE1), 
0.022 ± 0.002 mas yr _1 (0.108 ± 0.012 c, SW2-NE1) and 
-0.034 ±0.001 mas yr" 1 (0. 168 ± 0.005 c, J1-NE1). Positive 
velocities of SW1 and SW2 signify an advancing motion of 
both hotspots. The negative velocity of Jl indicates that 
the internal jet moves toward the terminal hotspot NE1. 
The separation of Jl at the first three epochs, 1995.266, 
1995.395, and 1996.375, shows a strong deviation from the 
general trend of the data points beyond 1997. During the 
fitting process, double uncertainties were assigned to these 
data points in order to reduce their weight and thus to avoid 
a bias induced by these large deviations. The separation of 
NE2-NE1 is about 1.2 mas, and the linear regression fit did 
not give a significant proper motion, because NE2 shows a 
significant transverse motion rather than a radial motion 

(Fig.EJ). 

Liu et al. (2000) reported a first determination of 
the NE1-SW1 angular separation rate of 0.058 ± 0.038 
mas yr -1 , derived from 8.4-GHz data with only six epochs 
between 1994 and 1997. Subsequently, Luo et al. (2007) 
improved this value to 0.031±0.006 mas/yr (5a), on the ba- 
sis of a longer time span (11 yr) and nine epochs. Compared 
to the previous 8.4-GHz work, the current proper motion 
measurements greatly improve the accuracy in the follow- 
ing ways: (1) The 15 GHz data provide a higher resolution 
(typically 0.8 x 0.6 mas 2 ), while the previous measurements 
at 8.4 GHz had a resolution (1.3 x 1.0 mas 2 ) that is com- 
parable with the separation between SW1 and SW2. At 
8.4-GHz the positions of SW1 and SW2 are affected by the 
flux density variation of the two components (mainly SW1), 
which is also true for the northeast lobe where the posi- 
tions of NE2 and Jl are affected by the flux density and 
structure variations of NE1. At 15 GHz, SW1 is clearly 
separated from SW2. (2) The present work covers a time 
span of 13.6 yr, four times longer than in the earliest work 
by Liu et al. (2000), and ~2.3 yr longer than that in Luo 
et al. (2007). (3) The data from 28 individual epochs give 
a better time sampling compared to the sparse sampling in 
previous works. (4) The 15-GHz VLBA data provide better 
and more uniform u,v coverage than the 8.4-GHz data. 

According to the new measurement of p = 0.027± 0.002 
mas yr" 1 for SW1-NE1 14ct detection), the kinematic age 
of OQ 208 is calculated as 255 ± 17 yr in the source rest 
frame. This age classifies OQ 208 as one of the youngest 



CSOs (iPolatidis k Conway (2003)1 |Gugliucciet al. (2005) 



Giroletti fc Polatidis (2009)[|An et al. (2012)p 



Assuming equal advancing velocities for NE1 and SW1, 
the internal jet Jl is found to move 3.5 times faster 
(vji = 0.23 c) than the terminal hotspots. A simi- 
lar phenomenon has also been obse rved in other CSOs 
(e.g., B0710 ±439 and B235 2+495: |Taylor et al. (2000)1 
J0132+5620 :IAn et al. (2012)[ ). While the radio lobes sweep 
the ambient ISM, a jet knot moves in an excavated channel 
with a relatively lower ISM density. 
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The kinematic parameters of the OQ 208 jet derived 
above can be used to calculate the jet flow properties using 
the following the equations: 



app 



r = 



/3sin(9 
1 - j3 cos ( 



f 



5 = 



1 



r(l - (3 cos 6) : 



(3) 



(4) 



(5) 



where, /3 app is the apparent speed of the jet component, 
j3 is the intrinsic jet speed (both in units of c), 6 is the 
inclination angle between the jet axis and the line of sight, 
r is the bulk Lorentz factor, and S is the Doppler boosting 
factor. Taking the apparent speed of the jet knot Jl (j3 apv = 
0.168 c) and the inclination angle (9 = 80.8°: Section 3.3) 
into account, we obtain f3 — 0.166 c, T — 1.014 and 8 = 
1.013. These calculations indicate a mildly relativistic flow 
in OQ 208. Relativistic beaming does not play a major 
role in determining the observed radiative and kinematic 
properties. 

3.5. Sideways motion and frustrated jets 

In addition to the advancing motion of the hotspots, 
OQ 208 shows evidence of sideways motion. Figure [U dis- 
plays the two-dimensional plots of the VLBI components' 
relative locations. The primary hotspot NE1 is used as the 
reference at the (0,0) position. The southwest hotspot SW1 
in general indicates an advancing motion along the NE1- 
SW1 line, although during the first half of the period 1996- 
2003 it shows a curved trajectory. Similar to SW1, the in- 
ternal jet Jl shows a general motion toward NE1, but it 
makes a loop-like path during this same initial period. This 
common feature of SW1 and Jl suggests that the reference 
NE1 may have followed a curved path between 1995 and 
2003. If SW1 is used as the reference (plot is not shown 
here), NE1 indeed follows a bending trajectory, first to the 
north and then to the east, in a mirror-symmetric pattern 
with SW1-NE1. Jl follows a straight trajectory along the 
connecting line J1-SW1. 

SW2 shows a complex and disordered motion pattern 
(Fig. |4j) , but generally it moves to the southwest along the 
connecting line NE 1-SW2 within the past fourteen years. 
Wang et al. (2003) measured a proper motion of SW2 (rel- 
ative to NE1) as 0.032±0.020 mas yr _1 based on five- 
epoch 5-GHz data. Their measurement agrees with ours 
(fj, = 0.022 ±0.002 mas yr _1 ) with a much higher accuracy. 

The hotspot NE2 shows an apparent motion to the 
southwest, in an opposite direction of the jet advance. As 
discussed in Section 3.1, NE2 is likely a hotspot generated 
by the deflected jet hitting the wall of the surrounding ISM. 
While the hotspot NE1 is advancing faster than NE2, NE2 
itself might be a stationary component. 

The wandering of the jet heads NE1 and SW1, as well 
as the disturbed lobe structure with deflected jet or double 
hotspots, provide the signature of an obstructed jet flow. 
According to the frustration model, the advance motion of 
the radio lobes is confined by the surrounding dense ISM 
due to the intrinsically low jet power, or the reduction or 
cessation of the jet power. A critical requirement for a CSO 



to evolve into medium-sized symmetric objects (MSO) is 
that the jet remains supersonic at the interface between 
the ISM and intergalactic medium (IGM). Through ana- 
lytic modeling of the expa nsio n of the hotspot and coco on, 
Kawakatu & Kino (2006) and Kawakatu et al. (2009) de- 
rived that the initial hotspot advance velocity should at 
least be 0.3c so that the CSO can evolve beyond the ISM- 
IGM boundary of the host galaxy. Jet flows with velocities 
below this threshold become subsonic before they reach 
the ISM-IGM boundary and have a distorted morphology. 

The present data show that the radio source of OQ 208 
is still growing, although the advancing velocity is low. The 
radio power of OQ 208 is Piaghz = 10 25 W Hz^ 1 with- 
out correction for the low-fre quency absorption, p lacing it 
in the low-jet-power regime (An & Baan (2012) I. A low- 
power CSO like OQ 208 can relatively easily develop hy- 
drodynamic surface instabilities in the jet, which makes 
the jet more likely to loss momentum. If a significant 
fraction of the jet momentum flux is dissipated, the jet can- 
not sustain the supersonic laminar flow and forms a stag- 
nating standing shock (i.e., the location of the hotspot) 
and becomes flaring and diffused beyond that point. Such 
a flared jet is found in the CSS quasar 3C 48, which has a 
compact an d bright hotspot about 3 parsec from the cen 
tral AGN 



(2 



Wilkin son et al. 
An et al. 



1991 
2010 



Worrall et al. (2004)| 
The decrease of jet 



Feng et al. (2005) 
power can happen in any stage of the radio source evolu 
tion when the nuclear activity is reduced or terminated, or 
the jet experiences signific ant loss of kinetic en ergy because 
of jet-ISM interactions ( |An fc Baan (2012)[ ). Frustrated 
sources may continue to grow but they do not have compact 
symmetric lobes. Eventually, frustrated CSOs and MSOs 
become radio relics at frequencies below a few hundred 
MHz. 



4. Conclusions and summary 

Multi-frequency multi-epoch radio images of OQ 208 with 
the highest angular resolution of ~ 0.5 mas present a typ- 
ical CSO morphology with two compact mini-lobes along 
a position angle about —126° and separated by ~10 pc. 
The flux density ratio of the structural components varies 
with the observing frequency. At 8.4- and 15- GHz images, 
each lobe is resolved into two subcomponents, among which 
the brightest components are identified as the primary 
hotspots. In the highest-resolution 15-GHz images, a knotty 
jet is detected between the NE and SW lobes. The core of 
OQ 208 cannot be securely identified from the present data. 

The two primary hotspots NE1 and SW1 show signif- 
icant flux density variations at 15 GHz of 62% and 19%, 
respectively. The 15-GHz variability of the hotspots in 
OQ 208 may result from the balance between feeding from 
the central engine and adiabatic expansion loss. 

The peak epoch of the NE1 flare in 1998.3 is ear- 
lier than that of the SW1 flare by an estimated 5.00 
yr, suggesting that NE1 is moving toward the observer and 
SW1 is receding. This light travel time difference between 
the two hotspots corresponds to a radial distance differ- 
ence of about 1.53 pc. Combining the projected and radial 
separation between NE1 and SW1, we estimate the inclina- 
tion angle between the radio jet and the line of sight to be 
~ 80.8°. This value is higher than the 70° opening angle 
generally assumed for the nuclear NLR, which suggest that 
the jet axis and the galaxy axis are not aligned. 
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Fig. 4. Relative right ascension and declination of the VLBI components. The primary hotspot NE1 serves as a reference 
position. 



Using the brightest hotspot NE1 as the reference, the 
relative proper motions were estimated for SW1, SW2, 
and Jl to be 0.027 mas yr" 1 , 0.022 mas yr" 1 , and -0.034 
mas yr _1 . The separation speed between SW1 and NE1 cor- 
responds to a hotspot advancing speed of 0.065 c, assuming 
symmetric advancing and receding hotspot motions. The 
proper motion of the jet component Jl relative to the sys- 
temic center is about —0.047 mas yr _1 and corresponds to 
a velocity of 0.230 c, making the jet mildly relativistic. The 
internal jet knot moves significantly faster than the termi- 
nal hotspots because it encounters a lower density in the 
excavated jet channel than the hotspots within the lobes. 
The angular separation rate leads to an estimate of the 
kinematic age of about 255 yr, suggesting that OQ 208 is 
one of the youngest CSOs known. 

The observed sideways motion of hotspots and dis- 
turbed lobe morphology are signatures of obstruction of 
the jet head by the surrounding ISM, although the overall 
radio source is still growing. During the early CSO stages 
of a radio source evolution, such disturbed lobe structures 



seem common (An et al. (2012) I. A young radio source may 
experience intermittent jet power or several failed starts be- 
fore entering into a continuous and steady growth phase. 
During each intermittent activity, the jets may impact dif- 
ferent regions of the ISM, failing to make a breakthrough. 
The notion of intermittent energy feeding is supported by 
the detection of multiple hotspots and the fading lobe ~30 
mas (~40 pc) to the southwest in a different position angle 
from the main jet body. Accordin g to the evolution mod el- 
ing of extragalactic radio sources (Kaiser & Best (2007) ), a 
critical requirement for a CSO to evolve beyond the ISM 
IGM boundary (typically 1-3 kpc) of the host galaxy is that 
the jet remains supersonic and maintains a laminar flow. 
Compared to high-power high-velocity sources, low-power 
low-velocity CSOs such as OQ 208 can easily develop tur- 
bulent jet flow after losing a significant amount of their 



momentum flux and kinetic energy during their interac- 
tions with the ambient ISM. Frustrated CSO radio sources 
with turbulent jet flow would not evolve into large-scale 
symmetric FR II- type radio galaxies. 
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Table 1. Model- fitting parameters 



Epoch 


S 


V 


Comp. 


R 


6 






P.A. 




(mJy) 


(GHz) 




(mas) 


(deg) 


(mas) 


(mas) 


(°) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(?) 


(8) 


(9) 


1995.266 


659.65(46.2) 


15.3 


NE1 


0.0 (0.000) 


-132.74 


0.365 


0.195 


0.28 




93.11 (6.5) 




NE2 


1.20(0.022) 


-26.32 


0.513 


0.513 


-35.54 




56.67 (4.0) 




Jl 


0.80(0.042) 


-148.46 


1.023 


1.023 


-112.62 




39.74 (2.8) 




SW1 


6.56(0.056) 


-125.95 


0.883 


0.883 


-149.04 




16.68 (1.2) 




SW2 


7.68(0.068) 


-116.73 


0.527 


0.527 


-107.10 


1995.395 


712.48(49.9) 


15.3 


NE1 


0.0 (0.000) 


43.77 


0.325 


0.155 


-4.56 




108.53(7.6) 




NE2 


1.15(0.035) 


-24.11 


0.549 


0.549 


-155.56 




42.50 (3.0) 




Jl 


0.77(0.044) 


-154.43 


0.400 


0.400 


-176.27 




39.84 (2.8) 




SW1 


6.58(0.070) 


-126.00 


0.855 


0.855 


141.34 




16.36 (1.2) 




SW2 


7.61(0.026) 


-116.87 


0.410 


0.410 


10.74 


1995.956 


723.35(50.6 


15.3 


NE1 


0.0 (0.000) 


-147.26 


0.484 


0.243 


-2.25 




75.48 (5.3) 




NE2 


1.24(0.037) 


-24.26 


0.594 


0.594 


-155.56 




21.71 (1.5) 




Jl 


1.45(0.062) 


-139.12 


0.631 


0.631 


-176.27 




36.42 (2.5) 




SW1 


6.54(0.074) 


-126.22 


1.011 


1.011 


-141.34 




15.22 (1.1) 




SW2 


7.66(0.084) 


-117.42 


0.454 


0.454 


-10.74 


1996.304 


780.52(54.6) 


15.3 


NE1 


0.0 (0.000) 


-132.59 


0.431 


0.199 


-9.72 




72.79 (5.1) 




NE2 


1.28(0.022) 


-24.40 


0.589 


0.589 


-155.56 




17.82 (1.2) 




Jl 


1.45(0.042) 


-139.21 


0.658 


0.658 


-176.27 




35.11 (2.5) 




SW1 


6.62(0.033) 


-126.29 


0.837 


0.837 


141.34 




15.50 (1.1) 




SW2 


7.67(0.045) 


-116.76 


0.537 


0.537 


10.74 


1996.375 


822.20(57.6) 


15.3 


NE1 


0.0 (0.000) 


33.49 


0.371 


0.179 


-11.31 




86.51 (6.1) 




NE2 


1.19(0.041) 


-25.63 


0.573 


0.573 


-59.04 




53.55 (3.7) 




Jl 


0.84(0.051) 


-150.52 


0.938 


0.938 


-114.44 




37.58 (2.6) 




SW1 


6.63(0.062) 


-126.42 


0.830 


0.830 


-176.19 




14.55 (1.0) 




SW2 


7.58(0.074) 


-116.57 


0.430 


0.430 


180.00 


1996.821 


767.11(53.7) 


15.3 


NE1 


0.0 (0.000) 


-5.09 


0.401 


0.167 


-13.44 




63.36 (4.4) 




NE2 


1.25(0.042) 


-24.97 


0.531 


0.531 


-26.57 




21.35 (1.4) 




Jl 


1.40(0.076) 


-137.89 


0.810 


0.810 


-162.90 




28.77 (1.8) 




SW1 


6.51(0.073) 


-126.84 


0.648 


0.648 


-120.96 




21.92 (1.5) 




SW2 


7.74(0.110) 


-117.57 


0.771 


0.771 


-122.01 


1997.656 


813.14(56.9) 


15.3 


NE1 


0.0 (0.000) 


4.54 


0.414 


0.187 


-13.52 




66.44 (4.6) 




NE2 


1.21(0.033) 


-26.95 


0.579 


0.579 


-155.56 




40.47 (2.8) 




Jl 


1.48(0.031) 


-140.71 


0.355 


0.355 


176.19 




33.54 (2.3) 




SW1 


6.59(0.040) 


-127.26 


0.706 


0.706 


141.34 




15.20 (1.0) 




SW2 


7.70(0.068) 


-117.68 


0.686 


0.686 


10.74 


1998.168 


864.91(60.5) 


15.3 


NE1 


0.0 (0.000) 


-128.99 


0.409 


0.182 


-17.91 




61.09 (4.3) 




NE2 


1.23(0.033) 


-27.39 


0.536 


0.536 


-141.34 




62.78 (4.4) 




Jl 


1.46(0.020) 


-142.44 


0.212 


0.212 


-118.61 




30.89 (2.2) 




SW1 


6.68(0.045) 


-127.20 


0.558 


0.558 


-156.80 




16.04 (1.1) 




SW2 


7.65(0.089) 


-117.18 


0.591 


0.591 


-145.31 


1998.179 


654.88(45.8) 


15.3 


NE1 


0.0 (0.000) 


129.29 


0.388 


0.153 


-16.93 




50.34 (3.5) 




NE2 


1.19(0.036) 


-28.80 


0.501 


0.501 


-71.57 




47.67 (3.3) 




Jl 


1.48(0.030) 


-142.30 


0.190 


0.190 


-156.04 




25.22 (1.8) 




SW1 


6.66(0.047) 


-126.95 


0.652 


0.652 


-140.71 




13.81 (1.0) 




SW2 


7.64(0.070) 


-118.03 


0.799 


0.799 


-118.61 


1998.234 


731.95(51.2) 


15.3 


NE1 


0.0 (0.000) 


-172.86 


0.441 


0.193 


-14.85 




50.07 (3.5) 




NE2 


1.23(0.035) 


-27.96 


0.570 


0.570 


-45.00 




57.65 (4.0) 




Jl 


1.46(0.025) 


-142.24 


0.250 


0.250 


-151.93 




25.95 (1.8) 




SW1 


6.65(0.042) 


-127.28 


0.594 


0.594 


-142.13 




16.05 (1.1) 




SW2 


7.72(0.074) 


-117.73 


0.791 


0.791 


-124.99 


1998.313 


909.06(63.6) 


15.3 


NE1 


0.0 (0.000) 


91.96 


0.419 


0.203 


-18.04 




67.76 (4.7) 




NE2 


1.11(0.027) 


-28.47 


0.549 


0.549 


-66.80 




72.93 (5.1) 




Jl 


1.47(0.017) 


-142.66 


0.240 


0.240 


-104.93 




40.75 (2.9) 




SW1 


6.73(0.038) 


-127.29 


0.775 


0.775 


-154.98 




13.36 (0.9) 




SW2 


7.76(0.085) 


-116.87 


0.605 


0.605 


-135.00 


1998.771 


853.14(59.7) 


15.3 


NE1 


0.0 (0.000) 


98.00 


0.412 


0.257 


-3.14 




52.21 (3.7) 




NE2 


1.20(0.027) 


-27.97 


0.331 


0.331 


-90.00 




97.50 (6.8) 




Jl 


1.42(0.016) 


— 144.00 


0.254 


0.254 


—61.39 




30.00 (2.1) 




SW1 


6.64(0.048) 


-127.83 


0.619 


0.619 


-142.43 




19.25 (1.3) 




SW2 


7.76(0.094) 


-117.94 


0.856 


0.856 


-105.52 


1998.828 


872.47(61.1) 


15.3 


NE1 


0.0 (0.000) 


164.27 


0.411 


0.194 


-17.25 




65.38 (4.6) 




NE2 


1.17(0.030) 


-29.16 


0.588 


0.588 


-129.81 




96.29 (6.7) 




Jl 


1.47(0.019) 


-143.83 


0.210 


0.210 


-119.75 




34.01 (2.4) 




SW1 


6.69(0.035) 


-127.61 


0.559 


0.559 


-153.44 




16.98 (1.2) 




SW2 


7.71(0.091) 


-117.78 


0.711 


0.711 


-144.46 
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1998.831 


861.22(60.3) 


15.3 


NE1 


0.0 ( 


0.000) 


90.43 


0.413 


0.189 


-19.03 




59.07 (4.1) 




NE2 


1.22 


;o.o4o) 


-29.05 


0.581 


0.581 


-155.56 




95.48 (6.7) 




Jl 


1.46 


;0.033) 


-143.81 


0.188 


0.188 


-176.27 




34.67 (2.4) 




SW1 


6.67 


;0.043) 


-127.61 


0.614 


0.614 


141.34 




16.32 (1.1) 




SW2 


7.68 


;o.o6i) 


-117.97 


0.611 


0.611 


10.74 


2000.991 


756.35(52.9) 


15.3 


NE1 


0.0 


0.000) 


-151.91 


0.393 


0.186 


-29.39 




47.06 (3.3) 




NE2 


1.20 


;0.032) 


-28.78 


0.541 


0.541 


-40.60 




94.08 (6.6) 




Jl 


1.42 


;o.oi6) 


-145.89 


0.178 


0.178 


173.66 




39.42 (2.8) 




SW1 


6.78 


;0.027) 


-127.37 


0.481 


0.481 


-165.96 




14.67 (1.0) 




SW2 


7.81 


;0.083) 


-118.44 


0.723 


0.723 


-112.62 


2001.010 


753.72(52.8) 


15.3 


NE1 


0.0 


0.000) 


124.68 


0.375 


0.195 


-31.21 




50.17 (3.5) 




NE2 


1.22 


;0.042) 


-29.31 


0.535 


0.535 


-107.10 




101.59(7.1) 




Jl 


1.40 


;0.028) 


-145.23 


0.213 


0.213 


-148.00 




37.17 (2.6) 




SW1 


6.74 


;0.039) 


-127.36 


0.438 


0.438 


-127.88 




17.37 (1.2) 




SW2 


7.77 


;0.099) 


-118.68 


0.824 


0.824 


-141.34 


2001.339 


758.44(53.1) 


15.3 


NE1 


0.0 


0.000) 


98.70 


0.344 


0.179 


-31.03 




57.32 (4.0) 




NE2 


1.09 


;0.045) 


-30.08 


0.538 


0.538 


-155.56 




106.86(7.5) 




Jl 


1.39 


;0.033) 


-145.15 


0.244 


0.244 


-176.27 




41.45 (2.9) 




SW1 


6.76 


;0.043) 


-127.09 


0.475 


0.475 


141.34 




16.03 (1.1) 




SW2 


7.71 


;o.m) 


-117.40 


0.684 


0.684 


10.74 


2001.785 


750.57(52.5) 


15.3 


NE1 


0.0 


0.000) 


-130.45 


0.357 


0.196 


-28.21 




57.76 (4.0) 




NE2 


1.13 


;0.029) 


-29.56 


0.547 


0.547 


-98.75 




110.58(7.7) 




Jl 


1.36 


;o.oi6) 


-144.15 


0.245 


0.245 


-163.30 




42.14 (2.9) 




SW1 


6.76 


;0.025) 


-126.95 


0.421 


0.421 


-156.80 




15.08 (1.1) 




SW2 


7.76 


;0.074) 


-117.31 


0.622 


0.622 


-142.43 


2003.313 


618.84(43.3) 


15.3 


NE1 


0.0 


0.000) 


-57.71 


0.395 


0.215 


-34.13 




54.58 (3.8) 




NE2 


1.14 


;0.053) 


-29.84 


0.616 


0.616 


-174.91 




125.72(8.8) 




Jl 


1.32 


;o.o3i) 


-144.22 


0.328 


0.328 


180.00 




42.00 (2.9) 




SW1 


6.78 


;0.046) 


-126.64 


0.524 


0.524 


153.92 




16.24 (1.1) 




SW2 


7.61 


;0.095) 


-117.74 


0.774 


0.774 


-167.73 


2005.473 


438.03(30.7) 


15.3 


NE1 


0.0 


0.000) 


-171.96 


0.383 


0.279 


-33.40 




59.50 (4.2) 




NE2 


1.15 


;0.042) 


-32.20 


0.616 


0.616 


-155.56 




136.50(9.6) 




Jl 


1.23 


;0.032) 


-142.98 


0.328 


0.328 


176.19 




40.62 (2.8) 




SW1 


6.85 


;0.042) 


-126.09 


0.524 


0.524 


141.34 




17.71 (1.2) 




SW2 


7.91 


;0.086) 


-117.47 


0.774 


0.774 


10.74 


2006.322 


358.16(25.1) 


15.3 


NE1 


0.0 


0.000) 


-103.38 


0.369 


0.232 


-26.88 




50.20 (3.5) 




NE2 


1.18 


;0.049) 


-33.44 


0.638 


0.638 


-111.80 




135.62(9.5) 




Jl 


1.18 


;0.035) 


-142.06 


0.341 


0.341 


168.69 




32.98 (2.3) 




SW1 


6.90 


;0.046) 


-126.14 


0.473 


0.473 


-146.31 




20.30 (1.4) 




SW2 


7.77 


;0.096) 


-117.99 


0.992 


0.992 


-106.70 


2007.493 


297.76(20.8) 


15.3 


NE1 


0.0 


0.000) 


144.95 


0.439 


0.342 


-12.76 




47.50 (3.3) 




NE2 


1.20 


;0.047) 


-35.96 


0.545 


0.545 


-172.56 




138.78(9.7) 




Jl 


1.14 


;0.033) 


-142.73 


0.337 


0.337 


-176.70 




35.89 (2.5) 




SW1 


6.91 


;0.048) 


-126.07 


0.590 


0.590 


-159.37 




17.50 (1.2) 




SW2 


7.92 


;0.078) 


-116.40 


0.738 


0.738 


-170.20 


2008.049 


290.45(20.3) 


15.3 


NE1 


0.0 


0.000) 


-172.66 


0.498 


0.405 


55.75 




42.41 (3.0) 




NE2 


1.24 


;0.068) 


-33.10 


0.680 


0.680 


-126.87 




117.24(8.2) 




Jl 


1.14 


;0.036) 


-143.25 


0.380 


0.380 


-168.69 




34.20 (2.4) 




SW1 


6.91 


;0.054) 


-125.97 


0.536 


0.536 


-158.20 




15.50 (1.1) 




SW2 


7.97 


;0.164) 


-117.49 


0.643 


0.643 


-126.87 


2008.131 


328.61(23.0) 


15.3 


NE1 


0.0 


0.000) 


114.90 


0.335 


0.335 


-161.09 




65.93 (4.6) 




NE2 


1.17 


;0.033) 


-33.18 


0.600 


0.600 


180.00 




160.05(11.2) 




Jl 


1.11 


;0.029) 


-142.00 


0.263 


0.263 


171.21 




39.06 (2.7) 




SW1 


6.96 


;0.034) 


-125.96 


0.561 


0.561 


180.00 




20.98 (1.5) 




SW2 


7.94 


;o.o5o) 


-117.06 


0.626 


0.626 


161.93 


2008.370 


267.87(18.8) 


15.3 


NE1 


0.0 


0.000) 


-161.26 


0.429 


0.355 


-8.18 




42.41 (3.0) 




NE2 


1.21 


;0.042) 


-35.83 


0.456 


0.456 


-172.56 




122.81(8.6) 




Jl 


1.11 


;o.oi6) 


-142.25 


0.277 


0.277 


-176.70 




31.06 (2.2) 




SW1 


6.96 


(0.04 ) 


-125.86 


0.506 


0.506 


-159.37 




17.25 (1.2) 




SW2 


7.94 


;o.m) 


-117.80 


0.764 


0.764 


-170.20 


2008.485 


264.86(18.5) 


15.3 


NE1 


0.0 


0.000) 


-123.69 


0.382 


0.382 


174.91 




49.90 (3.4) 




NE2 


1.16 


;0.043) 


-36.31 


0.553 


0.553 


-172.56 




135.60(9.4) 




Jl 


1.10 


;0.034) 


-142.38 


0.320 


0.320 


-176.70 




35.01 (2.4) 




SW1 


6.95 


;0.045) 


-125.97 


0.615 


0.615 


-159.37 




17.00 (1.2) 




SW2 


7.94 


[0.069) 


— 116.92 


0.699 


0.699 


-170.20 


2009.136 


216.27(15.3) 


15.3 


NE1 


0.0 


0.000) 


101.56 


0.403 


0.403 


69.90 




39.11 (3.0) 




NE2 


1.15 


;0.066) 


-30.72 


0.600 


0.600 


-155.56 




102.76(6.8) 




Jl 


1.08 


;0.035) 


-142.89 


0.320 


0.320 


176.19 




28.79 (2.0) 




SW1 


6.96 


;0.067) 


-125.39 


0.573 


0.573 


141.34 




14.65 (1.1) 




SW2 


7.93 


;0.123) 


-116.55 


0.713 


0.713 


10.74 
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2009.558 269.08(19.0) 15.3 NE1 0.0 (0.000) 118.00 0.463 0.36 48.97 

62.23 (4.4) NE2 1.13(0.041) -29.61 0.752 0.752 -47.49 

109.22(7.7) Jl 1.10(0.036) -142.84 0.326 0.326 33.74 

34.94 (2.4) SW1 6.91(0.042) -125.76 0.673 0.673 -86.25 

18.48 (1.3) SW2 7.89(0.058) -116.61 0.704 0.704 54.08 



